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Summary. We point out a possibility of the weak universality of spin-glass phase 
transitions in three dimensions. The Ising, the XY, and the Heisenberg models with 
± J bond distributions undergo finite-temperature phase transitions with a ratio of 
critical exponents 7/;^ = 2.4. Evaluated critical exponents agree with corresponding 
experimental results. The analyses are based upon the nonequilibrium relaxation 
from the paramagnetic state and the finite-time scaling. 



1 Introduction 

A spin glass (SG) phenomenon has been attracting great interests for thirty 
years both theoretically and experimentally [1] . It is also one of the toughest 
subject in a field of the computational physics. Though many theoretical in- 
vestigations have been done, it is still an unsettled issue whether or not a sim- 
ple random bond model explains the SG phase transition. Spins of many SG 
materials are well-approximated by the Heisenberg spins. However, numerical 
studies have suggested that there is no finite-temperature SG transition in the 
Heisenberg model [2,3]. This is quite different from the Ising SG case: many 
numerical investigations [4] consistently explain corresponding experimental 
results [5]. Kawamura [6] proposed the chirality mechanism in order to set- 
tle the discrepancy in the Heisenberg spin glass. The chiral- glass transition 
occurs without the spin-glass order. A small but finite random anisotropy in 
the real materials mixes the chirality and the spin. This effect induces the 
observed SG transition. The scenario sounds nice. However, Matsubara [7-9] 
recently recalculated the domain-wall excess energy and the spin-glass sus- 
ceptibility and suggested that the finite SG transition occurs. The methods 
are quite similar to the previous ones [2,3]. The final conclusion drastically 
changes by a subtle difference in the analyses of the obtained data. 

Computer simulations on SG models suffer from serious slow dynamics. 
The system sizes treated up to now are accordingly limited to very small ones, 
e.g., mostly a linear size L = 20 or less in three dimensions. The contradiction 
of the numerical results mentioned above possibly originates in the smallness 
of the system. This is a limit of the conventional equilibrium simulational 
study. One realizes the equilibrium state in finite size systems. Then, the 
thermodynamic properties are obtained by the finite-size scaling analysis. 
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We take an opposite approach to the thermodynamic Umit and reexamine 
the SG transition. First, wc take the infinite size limit by dealing with a very 
large system within a finite time range before the finite-size effect appears. 
Then, the finite-time scaling analysis is performed to obtain the thermody- 
namic properties. This method is known as the noncquilibrium relaxation 
(NER) method [10-12]. In our previous paper [13] it is clarified that the SG 
transition occurs at the same finite temperature as the chiral-glass transition 
occurs. Of course, the chirality freezes if the spin freezes. The estimated crit- 
ical exponent 7 is consistent with the corresponding experimental result [14] . 
The chirality mechanism is not necessary. In this article we expand the same 
analysis to the Ising SG model and the XY SG model in three dimensions. It 
is found that the Ising, the XY, and the Heisenberg models undergo finite- 
temperature SG transitions with the weak universality. 



The model we work with is the nearest-neighbor ±J random bond model. 
For each random bond configuration we prepare eight real replicas with dif- 
ferent paramagnetic spin configurations and with different random number 
sequences. Spins arc updated by the singlc-spin-fiip algorithm. We calculate 
at each time step t the spin-glass susceptibility Xsg through overlaps q'^^^ 
between two replicas (a and (3) out of eight: 



Indices n and v denote three components of spin degrees of freedom. The 

thermal average is replaced by the summation over the real replicas, whose 
number is denoted by m. It determines the statistics of the thermal average. 
Thus, a relaxation function Xsg(i) for one bond configuration is obtained. 
A configurational average is taken over Xsg{i) with different random bond 
configurations, different initial paramagnetic states, and different random 
number sequences. 

The spin-glass susceptibility is expected to diverge at the transition tem- 
perature (Tgg) as Xsg(i) ~ n/^'^, where z is the dynamic exponent. We obtain 
Tgg, 7, and zv by the finite-time scaling analysis. The NER of the Binder 
parameter g{t) is also calculated. Since the quantity is related to the fourth- 
order cumulant, quite a lot of samples are necessary to obtain meaningful 
data. Numbers of the bond configurations to obtain the results in this article 
is summarized in Table 1. The Binder parameter is expected to behave at Tgg 
as g{t) ~ i''/^, by which z is independently obtained. Then, v is estimated 
from a value of zv obtained by the scaling. The NER of the spin-glass suscep- 
tibility was originally applied to the Ising SG model by Huse [11] and that 
of the Binder parameter was applied by Blundell et al. [12]. 
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3 Results 



Numerical results suggesting the weak universality is shown in Fig. 1. The 
NER data of the spin-glass susceptibility exhibit an algebraic divergence with 
the exponent ^ j zv = 0.38. The NER function of the Ising model and that of 
the Heisenberg model are not distinguishable. That of the XY model differs by 
a factor three. An amplitude of the Xsg is dependent on the model. If we take 
into account a corrcction-to-scaling term with its exponent wjz = 0.5 [4], 
the relaxation functions fit to the expression At'^/^'^[1 — Bt~'^^^] from the 
first step. A deviation of the XY data from the diverging line after t = 10^ 
is possibly due to the finite-size effect. The NER of the Binder parameter 
shows an algebraic divergence with z = 6.2 for both models. Therefore, we 
obtain "f/f = 2.4 which is common to the Ising and the Heisenberg model. 
This suggests the weak universality. 
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Fig. 1. (Left) NER of the Xsg at Tsg. Data of the XY model are multiplied by 
three. The temperatures axe T/J= 1.17 for the Ising model, T/J = 0.46 for the 
XY model, and T/J = 0.21 for the Heisenberg model. The linear size, time steps, 
and numbers of bond configurations are summarized in Table 1. (Right) NER of 
the Binder parameter for the Ising model and the Heisenberg model 



Results on the Heisenberg model are briefly reported in [13]. We have 
obtained Tsg = 0.20(2), 7 = 1.9(5), and zv = 4.7(12). Since z = 6.2, v is 
estimated as v = 0.8(2). 

An example of a finite-time scaling analysis of the Ising model is shown in 
Fig. 2. This is an only preliminary result with a small system size {L — 19) 
and a short time range. The error bars are consequently large. Nevertheless, 
the obtained results are consistent with previous numerical investigations 
[4] and the corresponding experiment [5] as summarized in Table 2. The 
exponent estimated by the scaling {"i/zu = 0.39) is consistent with that of 
the raw data at Tgg as shown in Fig. 1 {^jzv = 0.38). 

Time evolution of the distribution function of the overlap P{q, t) is inves- 
tigated at T = Tgg = 1.17 (figures not shown). It exhibits a single Gaussian 
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Table 1. Numbers of bond configurations to obtain Xsg and g of Fig. 1 
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Fig. 2. (Left) NER of the Xsg of the Ising SG model at high temperatures. (Mid- 
dle) A finite-time scafing plot. We obtain r(T) for each temperature and 'y/zv 
so that the scaling is good. It is also possible for 7/01/ — 0.38 ~ 0.40. (Right) 
Estimate T^g supposing r(T) oc |T - T^gl''"'' 
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shape with a peak at g = before the size effect appears. As the time in- 
creases, the width of the Gaussian distribution grows in accordance with the 
divergence of the spin-glass susceptibility. It is possible to scale P{q, t) /t^/^^" 
plotted versus 

qt7/2z,. for various time steps from t=10to<=10^. That 
is, the NER function knows the critical phenomenon from the 10th step. 
The distribution changes its shape to having double peaks at ±qeq after the 
finite-size effect appears. 

It is also possible to investigate the low-tcmperature phase. The data sug- 
gest that the Xsg diverges algebraically with an exponent linearly dependent 
on the temperature. Details will be reported elsewhere [15]. 



4 Concluding Remcirks 

By applying the NER method it is clarified that the ± J spin glass models in 
three dimensions undergo finite-temperature SG transitions. There is possibly 
the weak universality with a ratio of the critical exponents j/i^ = 2.4. We 
compare our results with other numerical results and the experimental results 
in Table 2. They agree well within the numerical errors. Critical exponents 
7 and v in the Heisenberg model are almost one half of those in the Ising 
model. Reexamination of the experiment on the critical exponents is strongly 
encouraged. 
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The NER method is particularly effective in slow- dynamic systems. Ap- 
plications to various complex systems should be fruitful [16]. It potentially 
bring a breakthrough to a system stagnant due to the difSculty in computer 
simulations. 



Table 2. Estimates of Tsg, 7, v, /3, 7/1/, and z 
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Present work 


0.20(2) 


1.9(5) 0.8(2) 
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Matsuhara et al. [9] 


~ 0.18 
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Present work 


1.17(4) 


3.6(6) 1.5(3) 


~ 0.5 


~ 2.4 
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Mari-Campbell [4] 


1.195(15) 2.95(30) 1.35(10) 


~ 0.55 


~ 2.2 




Feo.sMno.sTiOs [5] 




4.0(3) ~ 1.7 


~ 0.54 


~ 2.4 
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